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We present the theoretical analysis of the steady state currents and density distributions of parti- 
, cles moving with Langevin dynamics, under the effects of an external potential displaced at constant 

rate. The Dynamic Density Functional (DDF) formalism is used to introduce the effects of the molec- 
ular interactions, from the equilibrium Helmholtz free energy density functional. We analyzed the 
generic form of the DDF for one-dimensional external potentials and the limits of strong and weak 
potential barriers. The ideal gas case is solved in a closed form for generic potentials and compared 
00 ' with the numerical results for hard-rods, with the exact equilibrium free energy. The results may 

CNJ , be of relevance for microfluidic devices, with colloidal particles moving along narrow channels, if 

external driving forces have to compete with the brownian fluctuations and the interaction forces of 
the particles. 
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INTRODUCTION 



Microfluidics, or the manipulation of liquids at the pico-liter scale, is an emerging technology with very interesting 
perspectives In particular, the construction of devices which might control colloidal fluxes at the level of single 

^ ■ particles, even under the effects of strong correlations at high density and in confined geometries would be of great 
interest for many applications, and it is already within the range of experimental works Q. A theoretical description, 
with quantitative prediction power, for these systems would allow the design of devices on a firmer basis, but such 
theoretical framework is well established only for the properties of systems at thermal equilibrium, mainly thought 
the applications of the density functional formalism (see e.g. In this paper we explore the use of a particular 

extension of the density functional methods to describe some systems out of equilibrium, like fluxes produced by the 
, time-dependent external forces in micro-pumps or turn-pikes which might control single colloidal particles in a flux. 
' The extension of the density functional formalism (DDF), from systems at thermodynamic equilibrium to the 
' dynamics of systems out of equilibrium has been a very attractive target for workers in classical and quantum fluids 
■ during the last two decades. The goal is certainly ambitious and the success has been limited to some particular 
' types of problems and levels of description. Thus, the modern theory of phase transitions dynamics Jlj is based on 
"j^ . quasi-local thermodynamic descriptions, which may be considered as an empirical extension of the squared gradients 
' density functional approximation to study the dynamics of a 'coarse grained' variable. The generic form for such 
descriptions, for a conserved 'coarse grained' density p{r, t) is 
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^ , where J-'[p] is the intrinsic the Helmholtz free energy functional, approximated by 

^ : T[p] ^ jdx (^f{p{x)) + \{^p{x)f^ , (2) 

X runs over all the space variables, f{p) gives the free energy per unit volume for bulk systems of homogeneous density 
p, and the coefficient h (which may be taken as a function of p{x) itself) gives a measure of the free energy associated 
to the inhomogeneity of the density. V{x) represents any external potential acting on the particles. 

The density functional approximation (jSJ has been extensively used to describe equilibrium systems with smooth 
spatial inhomogeneities, like the liquid- vapor interfaces J, 7], although it cannot describe systems with stronger 
modulations of the local density, like the layered structures of fluids near walls 0. The use of this T[p] for dynamics 
requires to specify two other empirical elements, the mobility r(p) which relates the density current to the local 
gradient of the thermodynamic chemical potential, and the random noise term S,{x,t) usually taken as a gaussian 
(uncorrelated) stochastic field which allows the systems go over thermodynamic barriers, e.g. in the process of 
nucleation of liquid drops from oversaturated vapors. Both the mobility as a function of the coarse grained density. 
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and the amplitude of the gaussian noise have to be taken as empirical mesoscopic parameters, which control the 
relaxation process of the coarse grained density, but without a direct connection to the molecular level. Therefore, and 
although equation ^ has been very useful to analyze different aspects of the phase transition dynamics , it cannot 
be used down to a molecular level of description, nor can it be extended to make use of the very accurate equilibrium 
density functional approximations ^ which go much beyond the simple functional form (0) in the description of 
equilibrium systems. Other approaches to define and use a DDF theory have used a more microscopic descriptions of 
the density distribution [3, llOl [ij but still keeping some coarse graining and a random noise term to overcome the 
energetic barriers of the coarse grained free-energy density functional. 

A different DDF approach, which a clear connection between the molecular hamiltonian and a dynamic density 
functional equation may be obtained |l2lll3| from the hypothesis of purely relaxative Langevin dynamics at molecular 
level; i.e. if the positions of the particles r,(t) (for i — 1,2, N) evolve in time according to the equation 



where i^(ri, rjv) is the potential energy which includes the interactions between the N particles in the system, and 
also any external potential acting on them. The molecular mobility To gives the (assumed constant) ratio between 
the mean velocity of the particles under a constant external force. When the theory is applied to colloidal particles 
moving in a fluid, the constant To may be obtained by Stokes formula in terms of their radius and the viscosity of the 
fluid matrix. The second term in the bracket in (O, represents a random force acting on particle i, with may be 
interpreted as the zero-average stochastic effect of the fluid matrix molecules on the colloidal particles. This is taken 
as gaussian uncorrelated function for each particle, but contrary to the random noise ^{x,t) in its amplitude is 
perfectly defined by Einstein relation, as twice the diffusion constant D = TokgT. The natural reduced units for the 
problem are kgT = 1//3 for the energy, molecular diameter a for distance, (3a'^/To for time, and hence To k^T/a for 
the velocity. 

Using the standard tools of stochastic analysis 0| and the equilibrium density functional formalism is possible to 
transform ^ into a Dynamic Density Functional (DDF) equation. 



within the following assumptions: 

i) The density p{x, t), following the dynamics of the DDF equation @, has to be interpreted as the average of the 
instantaneous density p{x,t) = X^^]^<5(x — i"i{t)), over an statistical ensemble of realizations for the random noise 
^(i) for the integration of the Langevin equation, from given initial conditions p{x,0) = Pi„i^{x). The temperature T 
is assumed constant and fixed by the bath. 

ii) The correlations between the colloidal particles in the dynamical system, at any instant t, are assumed to be 
identical to those in the equilibrium systems, which would be obtained by applied the appropriate external field to 
the system so that its equilibrium density distribution, p^c^ix) would be given precisely by the instantaneous value 
of p{x,t). These equilibrium correlations are those given by the density functional giving the intrinsic equilibrium 
Helmholtz free energy of the fluid, J-[p\. This equilibrium free energy functional, for classical fluids, has an explicit 
dependence from the ideal gas entropy contribution, 



and a contribution, AjF[p] = J-'[p\ — !Fid[p], from the particle interactions, which is known exactly only for a few 
systems; but approximations with increasing accuracy have been developed for others and give excellent results for 
the equilibrium properties in many cases. 

Beside having a perfectly defined mobility factor, as the density p{x, t) times the molecular mobility in the Langevin 
equation, instead of the generic function T{p(x,t)) in (^; the main difference between that semi-empirical equation 
and the DDF equation Q is the lack of any random noise term in later, which becomes a deterministic functional 
equation instead of the usual stochastic functional equation The averaging over the random force term in the 
original Langevin equation required to get p{x,t) with the interpretation i) above, gives the classical ideal gas 
contribution, .Fid[p] to J-^[p\ in the functional derivative of with the exact prefactor fc^T from the precise amplitude 
of the random force (t) ] Sind any extra random noise term in Q would produce the overestimation of the ideal gas 
entropy. The interactions between the particles, included in the derivative of Z// in Q give the contribution AJ-, to the 
intrinsic free energy functional J-[p] in Q, while the contribution to the potential energy from external fields appears 
explicitly as V{x, t). We may expect that the accuracy of the approach to include the effects of the interactions would 



dn{t) 
dt 



To ^tU{ri, ...,r7v) -l- 



(3) 




(4) 



^id[p] = kbT / dxp{x) {log[p{x)) - 1) 
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improve with the use of accurate forms of J-[p\, which include the non-focal functional dependence in a much better 
way than 

In particular, we may think of colloidal particles, interacting as hard-core molecules of diameter a and confined 
afong a hnear channel, as aheady explored experimentally i , which may be theoretically described as hard-rods 
in one dimension, for which the exact equilibrium free energy density functional is known ,,15.] . The solution of Q 
for hard- rods has already been studied for the relaxative dynamics under several static external potentials 
important conclusion of that work was that the use of the exact J^[p\ made unnecessary the inclusion of the random 
noise term in the DDF equation, since the exact free energy is always a convex functional, with only one minimum. 
The existence of activation barriers, between two different minima of the free energy, which should be overcome by 
the stochastic noise is an artifact of using an approximated form for J^[p]. Thus, very long relaxation times which in 
the solution of would be interpreted as typical of systems with a high free energy barrier between two minima, 
are recovered from the deterministic equation 10} without any free energy barrier, just from the particle correlations 
described through the exact free energy J-[p]- In this context, if we take colloidal particles of size cr « 10~^m, with 
mobility To given by Stokes law in water, at room temperature, the natural unit for velocity is ToksT/a w lO^^m/s, 
in a reasonable range for potential barrier shifts in microfluidic devices. 

In this paper we analyze the results of for a different class of problems, hard-rods particles in time-dependent 
external potentials V{x,t). In particular we start with systems in which the external potential moves with constant 
rate c with respect to the static framework in which equation (|2Jl applies. These problems would apply to colloidal 
particles moving along narrow channel under the force created by a moving electrical potential or optical twist; which 
may be used as a starting point for the design of more complex devices, like micropumps or turn-pikes. The plan of 
the paper is the following, in the next section we analyze the generic structure of the solutions of the DDF equation 
Q for the stationary states with density p{x, t) = p{x — ct) moving along the external potential V{x, t) = V{x — ct). 
In section 3 we get an explicit solution for the ideal gas case, when the interaction term Tip] may neglected; this 
solution helps to give a qualitative classification of the possible behaviour, in terms of the barriers of the external 
potential and the rate of displacement. In section 4 we explore, numerically the solutions of the DDF when the 
packing effects are important and described by the exact density functional at equilibrium for hard rods. We end 
with a general discussion on the observed qualitative behaviour and generic trend in the dependence with the physical 
parameters: the height and shape of the external potential barriers, the shift rate of these barriers with respect to 
the static framework for the relaxative dynamics and the effects of molecular packing. 



II. STATIONARY SOLUTIONS OF THE DDF 



The generic DDF equation Q may be interpreted as the exact continuity equation between the density of particles, 
p{x, t) and their current j{x, t), 

^ = -V,(.,0, (5) 
together with the approximated DDF closure for the current, 

j{x, t) = -To p{x, t) V + *0 ' 

in terms of the equilibrium free energy density functional. 

In the presence of a shifting external potential V{x, t) = V{x — ct) we look for stationary solutions of the form 
p{x, t) = p{x — ct), and j{x, t) = j{x — ct), so that the continuity equation Q implies that cVp(a;, t) = Vj(x, t), which 
for any problem in one dimension (ID) gives 

j{x,t)^cp{x,t)^A, (7) 

with an arbitrary integration constant A, which gives the difference between the actual current and that which would 
be produce by the rigid shift of the density profile at the same velocity c as the potential. This counter-drift term A 
cannot be determined from the DDF equation 10} , and it would be obtained only through the boundary conditions 
that we apply to the system. 

Using now equations © and ijT} we may obtain a closed functional equation for the stationary density p{x, t) = p(x), 
where we defined the variable x — x ~ ct, moving with the external potential V{x, t) = V{x); 
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In ID problems this expression may be formally integrated to give 



^ ^ + Vix) + — (x-Xo)-— —— = ^i, (9) 



Sp{x) To To Jx„ /5(£') 

where /i is the integration constant, linked to the arbitrary value of Xo in the lower limit of the integral. 

In the case of static external potentials, c = 0, the stationary density becomes the equilibrium density, and current 
j has to vanish everywhere, so that the counter-drift constant A in vanishes, and © becomes the usual Euler- 
Lagrange equation for the equilibrium density, p{x) — p{x), with the integration constant fi as the equilibrium 
chemical potential. Thus, in the dynamic case (c ^ 0) the integration constant /i in lO will still be used to control the 
total number of particles in the system, since for a given external potential V{x) and drifting rate c, that functional 
equation should have different solutions representing dynamic systems with different mean densities. The presence of 
a second integration constant A in is characteristics of dynamic systems, as it has to vanish in the limit c = 0. To 
determine the value of A we may consider a further symmetry assuming that the drifting potential V{x) and hence 
the stationary density p{x) are periodic functions with spacial period L (which imply a time period T = L/c). In 
that case, comparing the values of © for x and i + L we get 

To To Jx p{x') 



for any value of x, which implies 



and the stationary current is 

where the brackets represent the spacial average over a full period of the inverse local density. The last expression 
makes clear that for static periodic potential, c — 0, the only stationary solution of is the equilibrium state, with 
the current vanishing everywhere. 

The stationary density distributions and currents created by any periodic external potentials moving with a constant 
rate c may now be obtained from the solutions of ^ and , with the remaining integration constant p playing the 
same role as in the Euler-Lagrange equation, to control the mean density of the system. Before including an explicit 
form for the interaction between the particles, i.e. for the term AJ^[p] in free energy density functional, we may still 
get some generic analysis for particular limits. The case of a shifting external potential limited to a finite region of 
an infinite ID system may be obtained as the L — *■ oo limit of the periodic case. In that case, unless V{x) includes 
infinite barriers where p{x) = (as explored in section 2a below) we get A — cpo, determined by the density Po{p) of 
the uniform system at chemical potential /i, and equation may be written as: 

'^^^ 'n^)+^r d4i-^^=,- (13) 



5p{x) To y_oo V P(i') 

The current vanishes away from the region crossed by the external potential, and the integrated flow along the system 
becomes 

dxj{x) = c / dx{p{x) - po) = c{N - No), (14) 

i.e. c times the excess of particles generated by the shifting external potential, from the solution of (|13|) . This 
expression is just an integrated form of the continuity equation since the total transport of molecules is precisely the 
steady shift of the stationary density perturbation p{x) ~ po- 



A. Complete drift limit 



Let us consider a moving external potential which includes a periodic series of infinite barriers, where the density 
p{x) which solves © has to become zero. Then the average of the inverse density becomes infinite, whatever values 
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p{x) takes outside the barrier. The counter-drift term A vanishes and the current becomes equal to the density times 
the velocity of the external potential c. This is a regime of complete drift, where all the fluid particles trapped between 
two consecutive barriers of the external potential have to follow the movement of these barriers, so that j{x) = c p{x). 
Moreover, the functional equation ^ for the density distribution in that regime becomes equivalent to the equilibrium 
Euler-Lagrange equation 

^^ + V,ix)=p, (15) 

with a 'kinetic' external potential Vk{x) = V{x) + c/Fq x, i.e. with the true external potential plus a linear term 
proportional to c. In order to simplify the notation, we use a reduced drift rate with inverse length units, as c = (3c/To. 

The stationary density profile follows a variational principle equivalent to that of the equilibrium density functional 
formalism; i.e. p{x, t), described in terms of the drifted variable x = x — ct, is such that it minimizes the total 'kinetic' 
grand potential energy of the system 

n[p] = T[p] + J dxpix)iVk{x) - p). 

The constant slope c/Fo of the kinetic contribution to Vk{x) includes the extra force to keep the particles moving with 
mean velocity c, according with their Langevin dynamics © . The main advantage of this limit is that the solution of 
the DDF equation may be obtained with the well developed methods used for the equilibrium DFF, like conjugated 
gradients minimization of Q[p\. 

When the periodic shifting potential V{x) has high (but not infinite) barriers, we may expect that the stationary 
profiles p{x) would still be very close to the solutions of the equivalent equilibrium equation H15() . since everywhere 
outside the barriers p{x) 3> {p~^)~^ and the contribution from the integral in Q would be nearly constant (and hence 
transferable to p). However, the contribution of that integral over a full period has to exactly cancel the contribution 
of the linear term c i, to get a periodic equation, which is achieved by the integral giving a descending staircase 
shape, with smoothed steps at the barriers, and nearly planar plateau between them. Thus, with the appropriate 
choices of p in each case, the solutions of (|15|l over a single period, would be similar to those of the DDF equation 
everywhere except in the low density region at the high potential barriers, where the true stationary profile has 
to recover the periodicity. Only in the strict limit of infinite barriers the equivalent equation (|15|) may be regarded 
as exact, when used for the interval between point at consecutive barriers where p{x) = 0, as the role the integral is 
only to provide a constant shift of the chemical potential between each period. 

Moreover, even if the potential barriers of V{x) are apparently very high (on the fc^T scale) the quasi-complete 
drift regime described above may be frustrated when the drift rate c becomes too high, since the energy difference 
between the maximum and the minimum of Vk (i) diminishes with c, the effective barriers of the kinetic potential may 
become comparable with fc^T, and the counter-drift current grow. Thus, we expect a low velocity regime, in which 
the velocity of the particles v — j{x)/ p{x) is exactly c, and a saturation range in which further increase of c does not 
produce the increase of v. 



B. Weak external potentials 



We consider now the opposite limit when the external potential V{x) is very weak everywhere, so that it may be 
treated as a small perturbation over an otherwise uniform fluid. We may then expand the density around its mean 
value, p{x) — Po + Sp{x), with po — {p{x)), so that {l/p{x)) = p^^ — {5p{x)^) / pi + the current becomes now 

m^c(6pi£) + ^-^^+o\dp)y 

and its average is {j) = c {Sp{x)'^) / po + 0^{5p). Thus, a shifting weak external potential V{x — ct) = V{x) cannot 
produce the drift of the bulk fluid, the local current is proportional to c times the perturbation over the homogeneous 
fluid, with regions of positive and negative sign, and its average goes only as the density perturbation squared; i.e. 
the propagating waves in the density, p{x — ct) — p{x), produce a very small mean propagation of particles. 

So far in this subsection, we have only used the continuity equation and the requirement of having a periodic 
stationary density 1)10(1 of the DDF equation, to get the current in terms of the density perturbation. To get a 
relationship between the weak external potential, V{x), and the density perturbation Sp{x) we may use the expansion 
of the DDF equation Q around p{x) ~ po, as a generalization of the equilibrium linear response theory Q to the 
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dynamic case. In terms of the Fourier transform of the density perturbation 5p{q) and the external potential V{q), 
the linear term of ^ gives, 

[S{q,Po)] ^+ic/q 

where S{q, Po) is the equilibrium structure factor of the homogeneous system, which is proportional to the response 
function of the system to external field, as given by the second functional derivative of T[p\, by the standard DF 
analysis 0,^3- In the long wavelength limit it is given by the isothermal compressibility , S'(0,po) = hTpoX-riPo) 
while S{q,Po) = 1 in the short wave limit, q oo. The second term in the denominator of (|16|l comes from the drift 
and counter-drift terms in (Q, which give and imaginary part to the kinetic response function. 

If we consider a weak external potential with the form of an harmonic propagating wave, V(x — c t) — Vq cos{q x), 
the density would be p{x — c t) = po — pqCos{q x — 6), with amplitude and phase shift, 

P^q A a ^ -I ^il^Po) 

Pq = -Po Y/2 ^ 

\s{q,Po)]-^ + {c/qf] 

For low drift rates the density profile will be similar to that created by a static external potential, with the maximum 
density at the minima of V{x) and vice versa. When the external potential moves fast enough the effect of the 
imaginary part in the denominator of Hlt)|) puts the density perturbation out of phase with the external potential, 
and the maximum of p{x) moves toward the region with maximum external force, with 9 — > 7r/2, but the amplitude 
decreases as pg '--^ 1/c for large c. 

The mean velocity of the particles, in same reduced units as c, for any weak periodic potential may be written as 
a Fourier sum with the quadratic contributions of each Fourier component of the potential: 

ToPo ^[Siq,Po)]-^ + {c/qf 

The contribution of the short wavelength components, with c S{q, Po) ^ 9, is similar to the quasi-static limit 
(0 w 0), with V linear in c and quadratic with the static response S{q, po) V{q)] in the opposite limit, c S{q^ po) 3> g, 
the kinetic contribution dominates [9 ~ 7''/2) and the contribution to v is proportional to (V^(q)(7)^/c, i.e. decays 
with c, is independently of static structure factor and goes proportional to the mean squared force ((V^'(x))^), rather 
than to the mean squared potential. For each Fourier component of V{x), the threshold between the linear growth, 
"i; ~ c, for low c and the decay, v ~ 1/c, for large c, would give the optimum shift rate c^^^ q/S{q,po) to get the 
maximum mean velocity. In the ideal gas limit, S{q, po) = 1, this optimum value is c = g, becoming very small in the 
long wavelength limit. For correlated systems, there would be a non monotonic dependence of the maximum v with 
c, depending on q due to the static structure factor. 



III. SOLUTION OF THE DDF EQUATION FOR THE IDEAL GAS 

In the previous section we have obtained some generic limiting behaviour for the solutions of the DDF equation. 
We turn now to explore the general solution of 1^ , first in absence of particle interactions, and in the following section 
including the effects of molecular packing as one dimension hard rods. 

The ideal gas case corresponds to take AjF[p] — and it should be the exact low density limit of any fluid system. 
Moreover, the assumptions to get Q from Q become irrelevant, and the DDF equation may be consider as an exact 
Fokker-Planck equation for the density distribution. The equation JSJ becomes 

Vp{x)+p{x){pVV(x)+c)^^. (18) 

When A is zero (i.e. in the static limit, c = 0, or in the complete drift limit described above), the homogeneous 
linear differential equation is directly solved to give p{x) — A exp{— j3Vk{x)) in terms of the Boltzmann factor of the 
kinetic potential (3Vk{x) = l3V(x) + cx. The arbitrary constant A is the activity, playing the same role as the chemical 
potential /i in lO to control the total number of particles. 

The general solution of the inhomogeneous equation with A 7^ is obtained, by the usual Green function method. 
For periodic external potentials, with period L, it has the generic form, 

p{x) = A e-^^^(^) f 1 + ~ ^ r dy e^^'-'^y^] , (19) 



dy el^^'-iy) Jo 
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which may be evaluated in terms of single integrals of the inverse Boltzmann factor of Vk (x) . The counter-drift term 
A, is To A/ (3 times the prefactor of the integral, inside the bracket. 

The average velocity of the particles, defined as the ratio between the mean current {j{x)) and the mean density 
{p{x)), and in the same reduced units as c is 

- ^^'^'^^ -<^^, (20) 



with the dimensionless auxiliary function 




dy (^e 



t3{V{v)~V{x)) _ -^^ gCy 



dye^^^y^+'y] (21) 

/ cL 

It is instructive to analyze the structure of this equation for potential barrier V{x) which is zero outside an interval, 
taken as < x < a, much shorter than the periodic length of the system L. When the maximum height of the 
potential barrier, pV{x), is well above the maximum kinetic term cL, the results of H19I20|I is always within the 
quasi-complete drift regime described above, with A ^ c (p) and w « c. For larger drifting rates, when cL becomes 
much larger than the maximum of pV{x), the main contribution to ^(c) comes always from the first term in the right 
hand side of H21|) . with a delicate balance between the fast growing exponentials in the integrals. 

The simplest case of a square barrier with V{x) = Vo in an interval of width a and V{x) = outside, may be solved 
analytically to give ^(c) and v. For large drifting rates, c mm[a, (L — a)] ^ 1, we get 

- - ' (22) 



cL + el^^" - 1 ■ 



Therefore, a square potential barrier of height Vo would be enough to produce the complete shift of the ideal particles 
up to (reduced) drifting rates c well below {e^^° — 1)/L. For c beyond this limit the mean velocity would saturate 
to a constant value v « e^^°/L for (3Vo ^ 1 and v ~ PVo/L for [3Vo ^1. In Figure 1 we compare the values 
of V for square barriers, parabolic barriers, V(x) = Voil — {2x/w)'^) for —w/2 < x < w/2, and gaussian barriers 
V{x) — Vo&x'p\—x^ I (2c?")\. At low c the three barriers give similar results, with v proportional to c both for high 
barriers (Fig. 1(a)) and low barriers (Fig. 1(b)). The fortes case corresponds to the complete drift regime u = c, 
predicted in section 2(a) for large enough values of fSVo', while the second corresponds to the linear response result 
V ^ c {JJVoY 1^^ predicted in section 2(b) for low c. At higher values of the shifting rate c, the results of the square 
barrier differ qualitatively from those of the parabolic and gaussian barriers. The two later reach a maximum value 
of the mean velocity, v„iax = v{cjnax), and then decay as v ^ 1/c, while the square barrier produce a monotonous 
increase of 5 up to a constant value for large c. The comparison of the results with different barrier shapes and 
parameters shows that the most relevant aspect of V{x) is the maximum slope, rather than the barrier width or 
the barrier height. Thus, reducing the width of a parabolic barrier, while keeping unchanged its maximum value 
Vo, increases rather than decreases the maximum value of the mean velocity, since it increases the maximum force 
acting on the particles. In this sense the square barrier model should be consider unphysical, since the infinite force 
acting at the barrier borders creates the saturation of v at large c. Choosing the parameters of the parabolic and 
gaussian barriers to have similar maximum slopes, produces very similar results for w(c), as shown in Figure 1 for the 
empirical choice oiw = 5.5a and the same barriers height, Vo- This robustness of the results to the particular form of 
V{x), as far as it does nor include unrealistic discontinuities, might help to eventual comparisons of the theory with 
experimental results. 

The weak potential limit H17|) . integrated over the Fourier components of the different barriers V{x), may be 
compared with the exact results obtained from H20I21() . The qualitative difference between the sharp square barrier 
and the gaussian or parabolic models is a consequence of the slow, V{q) ~ 1/q, decay of the external potential Fourier 
components for the discontinuous square barrier; while the first derivative discontinuity (i.e. V{q) ~ 1/?^) for the 
parabolic barrier is already weak enough not to change the qualitative dependence of v{c) with respect to the smooth 
gaussian barrier. For the later, and assuming that the distance L between barriers is large compared with the barrier 
width a, we may perform the sum over Fourier components in l|17|l as an integral to get an analytic form for the mean 
velocity within the linear response approximation: 

^(c) = ac{^\pii — nac e"^ Erfc[ac]), 
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in terms of the complementary error function which leads to 

I 

L \2ac 4(ac)2 

for large c and 



^(c) = ? TT-, 



v{c) = ^ (ac - ^[a-cf + 3(ac)3 + ....) , 

low c; with a maximum Wmaa; = 0.4346(/3Vo)^/i at Cmaj; = 0.82132/a. 

The results in the inset of Figure 1(b), shown that for f3Vo = 1 the linear response approximation for v{c) is 
very similar to the exact result. In general, we may expect such good agreement for gaussian barriers up to barrier 
height (3Vo ~ (a/L)~^/^, beyond this limit the linear response for low drifting rates would predict unphysical values 
of v{c) > c, beyond the complete drift limit; which would indicate the breakdown of such linear response analysis 
since would predict negative densities at the potential barrier. However, for large drifting rates, when v ~ 1/c, the 
predictions of the linear response approximation are quite good even for very large values of (3Vo, as shown in the 
inset of Figure 1(a). The optimum value c = Cmax to get the maximum v would be independent of pVo in the 
linear response approach, while the exact result gives Cmax growing roughly linearly with /3Vo, as the indicated by 
the crossover between the regimes v ^ c and v ~ (pVo)"^ / {Lac) , for low and high values of c, respectively. However, 
the maximum value of v from the exact result remains rather close to the predicted by the linear response, even at 
barriers as high as pVo = 20, despite the different predictions for the value Cmax at which that Vmax is achieved. This 
behaviour is clearly shown in Figure 2, with logarithmic scales for both pVo, Vmax, Cmax and the slope of v(c) at c = 0, 
and it may be explained as the results of cancellations of two different effects: a tendency of the exact v{c) to grow 
faster with pVo than the prediction of the linear response, for low barriers, and the tendency of Cmax to grow linearly 
with pVo at high barriers, leading to Vmax ~ Cmax ~ pVo/VaL for very large (3Vo. Similar results, with different 
analytic forms for v(c) may be obtained for the parabolic barriers or other parametric forms of the external potential. 

In Figure 3 we present the density profiles and the steady currents for a smooth gaussian barrier of width a — 0.707cr, 
space period L = 30a, and two set of the barrier heights and the shifting rates. For /3Vo — 1 the results of the linear 
response theory gives a reasonably accurate description of the exact density profiles for the slowest shifting rate, 
c = Icr^^ in Fig. 3(a) and a very accurate results at c = lOcr^^ in Fig. 3(b). The shape of p{x), relative to 
the potential barrier (shown by the dot-filled gaussians) reflects the difference between the quasi-static behavior for 
c = 1<7~^, and the strong effects of the kinetic effective potential Vk(x) = V(x) + cx a.t c = lOti"^, which pushes the 
maximum of p{x) toward the maximum of V{x). The current, which may be represented by the same line with the 
appropriate changes in the vertical axis to get lO, has a positive contribution in front of the barrier partially cancel 
by the negative contribution across the barrier. 

For higher barriers , f3Vo — 5, and low drift rate Figures 3(c), the system becomes close to the complete drift 
regime (A — > 0), with a strong accumulation of particles in front of the advancing barrier, which decays exponentially 
p{x) exp{—cx), so that for long distances L between the periodic barriers the density at the receding side is 
extremely low compared with the density at the advancing side. The linear response results are quite out of the 
mark, predicting unphysical negative density at the barrier. At higher values of c Figures 3(d) the accumulation of 
particles at the barrier front is so large that the barrier becomes ineffective to keep the complete drift regime, the 
local chemical potential, p,{p) = f3~^log{p) becomes comparable to Vo and the particles jump over the barrier. The 
final effect is that the counter-drift A grows and the total effective potential in ^ becomes much weaker that the 
external potential V{x). In this way, the current decreases and it becomes closer to the linear response result, with 
the limit of no density perturbation or steady current for very rapid displacements of the barrier ( i.e. c — > oo), for 
any smooth shape of the barrier. 



IV. STEADY STATE FOR HARD-ROD PARTICLES 



In the solution of the DDF equations for the ideal gas case the mean density of the system is irrelevant for the mean 
velocity v, as both the density p{x) and the current j{x) have the common factor A to include the constrain of fixed 
number of particles. In this section we explore the effects of molecular packing and correlations on the stationary 
currents created by the shifting external potentials V{x — c t) = V{x). For particles moving along a linear channel, 
the choice of hard-rods (with length a) to model the interactions allows the use of the exact DF for the equilibrium 
free energy [Tsj l . The functional derivative of the equilibrium free energy in Q is 



(23) 
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with the local packing fraction ri(x) defined as 



dx p(x'). 



(24) 



The equation© for the stationary density becomes now an integral equation, to be solved with periodic boundary 
conditions and the the self-consistency requirement The total number of particles in a period = L{p{x)), 

should be controlled by the chemical potential fi, with non-trivial dependence of p{x), A and j{x). Alternatively, we 
may solve dynamic equation for p(x,t), from a reasonable guess for the initial density distribution at t = 0, and 
integrate in time until the stationary limit p{x, t) = p(x — c t) is achieved. In that case, the number of particles is 
set with the initial guess (since equation (@J keeps constant N). We have checked that both procedures lead to the 
same stationary results, and the relative computational efficiency would depend on the characteristics of the external 
potential, and the choice of the initial profile. 

For weak perturbations (weak external potentials or very large drifting rates) the linear response analysis 1)161) . has 
to be applied with the hard-rods structure factor, 



with the packing fraction rj = apo- Near the complete packing limit t] ~^ 1, the effect of a shifting oscillatory external 
potential V{x,t) — VqCos{q{x — ci)), would depend strongly on the wave vector qa; as shown in Figure 4(a) for 
constant value of Vq (i.e. constant amplitude for the oscillation of the external potential), the optimal choice of q to 
get maximum mean velocity v, at low values of c, corresponds to q at the first peak of the structure factor, just below 
qa = 27r, with a large enhancement over the velocity obtained in an ideal gas with the same external potential. For 
larger values of c the maximum v is obtained for q in the second or in the successive peaks of S{q,po)- The results 
in Figure 4(b) for constant qVq (i.e. constant amplitude of the oscillating force) indicate that the optimum choice in 
that case is always at the first peak of the structure factor, but the enhancement over the ideal gas is strongly reduced 
as c increases. 

For well separated periodic barriers, like gaussians of width a much larger than their period L, the main effects 
of the hard-rod interactions are to decrease the compressibility and to form layered structures on both sides of the 
barrier, as usual in the equilibrium density distributions of a hard-core fluid near a wall. In figure 5 we present the 
results for a dense system, with mean density N/L = 0.8a~^, under the action of high gaussian barriers {PVo = 10.6 
and width a = 0.707(t, separated by a distance L — SOa. At low shifting rate (c = 0.2a~^ in Fig. 5(a) ) the system 
is close to the complete drift limit, with very low density at the potential barrier. On both sides of the barrier here 
are strong oscillatory structures, around a nearly constant density. The amplitude of the oscillations is larger at the 
advancing side than at the receding side of the potential, but this may be understood in terms of the equilibrium 
density profiles for the effective kinetic potential Vk{x) = V{x) + (3~^cx, which makes sharper the effective wall at 
one side and smoother at the opposite side. 

Increasing c, in Figures 5(b-c), produces an apparently weak effect on the density at the outside the barrier. For 
comparison we include the ideal gas results, for the same mean density and barrier (3Vo = 5, chosen to be equivalent to 
the effective barrier in the hard rod, with the recipe describe belong , which have high peak densities at the advancing 
front with maximum values of Pmax = 3.3cr~^ for c = la~^ and/Jmax — 7.5a~^ for c = 1.7a~^, with exponential 
decays toward a much lower on the receding edge of the next barrier. The low compressibility of the dense hard-rod 
fluid strongly reduces such accumulation of particles in front of the barrier and their sensibility to the value of c; the 
main change between the results for hard- rods in Figures 5(a) and 5(b) is in the increase of the minimum density , 
near the barrier maximum, which indicates that the local chemical potential may get easily over Vo with very weak 
changes in the local density, as the compressibility of the system {(3 p)'^^ dp / dp = (1 — 77)^ produces strong changes 
in p with small changes in ?/, and the particles start jumping over the barrier without building up the high local 
densities of the ideal gas. Only for the largest value of the shifting rate, c = lOfi"^ in Fig. 5(d), we start to see a clear 
accumulation of particles at the advancing front, but still being a small fraction of the total number of particles in the 
system. The oscillatory structures become damped, with respect to their shapes for lower c and this effect is stronger 
on the advancing than on the receding front, as a signature that the effect of term proportional to A in the effective 
potential of the Euler-Lagrange equation ||^ is getting over linear kinetic contribution ex. The linear response results, 
presented as the dotted lines in Figures 5, are quite different from the full DDF calculations, as could be expected for 
such a high barrier, except for the largest value of c = 10cr~^. However, even in that case the linear response fails to 
reproduce the layering structure. 

Finally, we present in Figure 6, the mean velocities produced by gaussian barriers, as functions of the shifting rate, 
in system of different mean density. The results for the ideal gas with a barrier height of l3Vo = 5 are compared 
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with those for mean densities of Po = N/L = 0.2, 0.4, 0.6 and 0.8cr~^, with barrier heights adjusted to give the same 
effective barrier, Vo — A^(po), in terms of the excess chemical potential over the ideal gas, at each mean density. This 
semi-empirical recipe is effective to give similar values of v for different po, while keeping the same absolute barrier 
height produces a rapid decrease of v with the mean density. In all the cases the qualitative shape of v{c) is similar 
to that of the ideal gas, with a linear growth, w « c in the quasi-complete drift regime, followed by a maximum mean 
velocity and a decrease w ~ 1/c at large shifting rates. The particular values of Vmax are tied to our empirical choice 
of the values of (3Vo at each mean density, but the optimal rate Cmax at which it is achieved, and the generic shape of 
the curve are roughly independent of the barrier height. In this respect, the main effect of the hard-core interactions 
seems to be the broadening of the maxima which becomes much less sensitive to the value of c than for ideal gas 
systems. For very large shifting rates the linear response analysis becomes valid again, as for ideal gas systems, and it 
becomes independent of S{q, p) and hence equal to that for the ideal gas. This may be interpreted, in physical terms, 
as the result of a very rapid series of of pushing and pulling forces, acting on the particles as the potential barriers go 
over them, with a time period much shorter than the typical collision time between particles; so that the later become 
irrelevant. 



V. CONCLUSIONS 

We have applied the dynamic density functional (DDF) theory to study the relaxative dynamic of particles moving 
along one-dimensional channels, under the action of external potential shifting at constant rate. Our results could be 
of interest for systems of colloidal particles, with typical sizes of cr ~ 10~^m, following brownian motion on water at 
room temperature. In that case the relevant values of the shifting rates of the external potential would be around 
10~'^m/s, which could be achieved by the the dielectric force created on the colloidal particles by a moving laser beam 
or the force created on surface charged colloidal particles by a moving electric potential. Although strongly simplified 
with respect to the complexities of real micro-fluidic devices, we may hope that the model and its theoretical treatment 
contains the essential features of real systems. In principle the model could be studied directly from the stochastic 
equation for the brownian motion of the particles, without relaying on the equilibrium density functional formalism 
and its approximate application to the description of systems out of equilibrium. For the ideal gas case the results 
should be identical, but the DDF theory allows a direct extension to interacting systems, which we have considered 
here as pure repulsive interactions treating the colloidal particles in a narrow channel as hard rods particles in one 
dimension. 

Under the effects of a periodic external potential shifted at constant rate c, the systems would achieve a steady 
state, with density distributions and local currents following the shift of the external potential. These are the only 
problems analyzed in this paper, leaving for future works the analysis of the transient evolutions to achieve such steady 
state from arbitrary initial conditions. The steady states are particularly appropriate to be studied within the DDF, 
since the structure of the relevant equations becomes similar to that of the usual density functional formalism for the 
equilibrium properties. Thus, the generic equation Q is equivalent to an Euler-Lagrange equation for the equilibrium 
density distribution but including two extra terms as external potentials. The first one, which we have called the 
kinetic contribution, is just a linear potential, proportional to the shifting rate c which breaks the symmetry between 
the advancing and the receding front of a shifting potential barrier. The second contribution to the effective external 
potential is a self- consistent term which depends on the integral of the inverse density distribution over the system, 
and it is proportional to the second integration constant. A, which (together with the chemical potential p) is required 
to get the general solution of the integro-differential Eulcr-Lagrange equation for the problem. For external potentials 
made of periodic barriers, we have given an intuitive interpretation to A as the counter- drift current, associated 
to those particles which jump over the barriers and so reduce the mean current in the system, with respect to the 
complete- drift regime in which the full density is shifted at the same rate as the barriers, and the density distribution 
may be obtained as the equilibrium profiles for the fluid in presence of the external plus kinetic potentials. This 
regime would appear, and hence the integration constant A would vanish, whenever the potential barriers are large 
enough, compared with the kinetic effective potential. In the opposite limit, for very weak shifting external potentials, 
we have extended to the DDF the usual linear response analysis of the equilibrium density functional formalism. The 
range of application for such analysis is in fact enlarged with respect to the equilibrium case, since strong external 
potentials moving at very large rate would produce much weaker effect than the same external potentials at rest. In 
that limit, opposite to the complete-drift, the self-consistent contribution to the effective potentials cancels out most 
of the linear kinetic term and a good part of the external potential, so that the density distribution is much smoother 
as a function of the position that the static equilibrium density profile. Only for discontinuous external potentials, 
like square barriers, there would be a remnant effect of the shifting external potential in the limit of c ^ oo, and it 
should be considered as an artifact of the infinite forces acting on the particles at the jumps of V{x). 

Although we have used periodic potential barriers in our analytical and numerical solution of DDF problems, we 
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have also considered the case of single potential barriers as the L ^ oo limit of the periodic potentials. In that 
case, the integration constant A goes to the shifting rate c times the bulk density of the fluid, Po far away from the 
moving external potential. This is the requirement to get null current away from the external perturbation as the only 
possible solution for the stationary density distribution in an infinite one-dimensional system with a shifting potential 
acting on a finite region. Any solution of the Euler-Lagrange equation with A 7^ cpo would imply the accumulation of 
particles at the advancing front of the potential barrier and the depletion of the density behind it, so that the system 
would not be at any stationary state. In any realistic problem, like possible devices to control the flow of colloidal 
particles along channels, the finite length of the channels would probably produce typical density distributions and 
currents closer to those analyzed here for periodic systems than to those predicted for infinite long channels. 

For the periodic systems studied here, we have found that the mean velocity of the particles may be maximized 
by the appropriate choice of the shifting rate. The hard core interactions among the particles reduces the optimal 
value of Vmax, but makes this mean velocity a broader function of the shifting rate. At high packing fractions, Vmax 
has a strong and non-monotonic dependence on the period of the external potential, as a results of the correlations 
between particles. Such effect could be used to design more effective micro-pumps or turn-pikes, in which the flow of 
colloidal particles would be controlled by the oscillatory external potential. The particular models studied here, ideal 
gas and hard-rods, under the action of gaussian potential barriers or harmonic periodic potentials, have been chosen 
to explore the main regimes and generic properties of stationary flows in one-dimension. There are obviously many 
open questions to answer and different models to be explored. The application of the density functional techniques to 
these systems is still at its beginnings, and it would take time to ascertain its full capabilities. However, the generic 
theoretical requirements: the absence of hydrodynamic effects and the similarity between the particle correlations in 
the dynamic and in the equivalent equilibrium systems, seem likely to be reasonable assumptions in systems of colloidal 
particles moving at rates of the order of microns per second, under the effects of controlled external potentials. The 
application of the DDF techniques to systems in higher dimensions and to non-stationary problems would hopefully 
benefit from what we have learned with our simpler models and problems. The design of practical micro-fluiditic 
devices could also be supported by the generic principles developed with these simple models, and the experimental 
realization of these systems appears to be within the range of the present technologies. 
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FIG. 1: Mean velocity of the particles as a function of the barrier drift rate, for square barrier (solid line) of width a = 2cr; 
parabolic barrier (dotted line) of width w = 5.5a and gaussian barrier (dashed) of width a = la. The barrier heights are 
/3Vo = 7 in figure (a) and pVo = 1 in figure (b).The insets compare, for the gaussian barriers, the exact results (solid lines) 
with the linear response theory (dashed lines). 



FIG. 2: Maximum mean velocity value Vmax, drift rate, Cmax, at with this maximum mean velocity value is achieved, and the 
derivative of mean velocity at zero drift rate, v'(ac = 0) , for a gaussian barrier as functions of the barrier height. The solid lines 
give the exact values and the dashed lines give linear response approximation, which are similar to the exact for low barriers. 



FIG. 4: Mean velocity of the particles under a shifting oscillatory external potential as a function of its wave-number q, for 
constant value of the barrier amplitude, /3Vq = 1, in (a), and constant force amplitude qa/3Vq = 1 in (b). The full line gives 
the result for hard-rods at po = O.Scr"^ with c = 2a~ , and the dashed lines the same for c = Wa~^. The dotted lines gives 
the ideal gas (i.e. poO" = 0) under the same drifting potentials for c = 2a~^ and the dash-dotted line for c = lOu"^. All the 
results are calculated with the linear response theory. 



FIG. 5: Density distributions (solid lines referred to the left axis) and currents (the same lines but referred to the right axis), 
for hard rods of mean density po = 0.8cr~^ under the infiuence of an external gaussian potential of height ffVo = 10.6 and width 
a = 0.707(7. The density distribution predicted by the linear response theory (dotted lines) and the ideal gas for an external 
gaussian potential of height pVo ~ 5 and width a = 0.707(7 (dashed lines) are also presented for comparison. The drift rates 
are c = 0.2a~^ in (a), c = l.Oa"^ in (b), c = 1.7a~^ in (c), and c = lOcr"^ in (d). 



FIG. 6: Mean velocity of hard-rod particles, as function of the drift rate, c, produced by a gaussian barrier of high 0Vo = 5 in 
an ideal gas system, and a gaussian barriers of height fiVo ~ 5.47, 6.18, 7.41, and 10.6 for systems with mean density PoC = 0.2, 
0.4, 0.6, 0.8 respectively. The width of all gaussian is a = 0.707it. 



FIG. 3: Exact density distributions for an ideal gas system (solid lines referred to the left axis) and currents (the same lines 
but referred to the right axis), for gaussian barriers of width a — 0.707(7 and heights f3Vo = 1(a) and (b) and /3Vo = 7(c) and 
(d) . The drift rates are c = 0.5a^^ for (a) and (c), and c = Wa~^ for (b) and (d). The dot-filled lines represent the external 
potential barrier (not at scale) just to show the relative position of the density structures to the barrier. The dashed lines 
(referred to the left axis) give the density distributions predicted by the linear response theory. 



